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Abstract. We present a summary of our recent rigorous derivation of the cele- 
brated Ginzburg-Landau (GL) theory, starting from the microscopic Bardeen- 
Cooper-Schrieffer (BCS) model. Close to the critical temperature, GL arises 
as an effective theory on the macroscopic scale. The relevant scaling limit 
is scmiclassical in nature, and semiclassical analysis, with minimal regularity 
assumptions, plays an important part in our proof. 



The purpose of this paper is to describe how the celebrated Ginzburg-Landau 
(GL) model of superconductivity 1 arises as an asymptotic limit of the microscopic 
Bardeen-Cooper-Schrieffer (BCS) model [2]. The relevant asymptotic limit may be 
seen as a semiclassical limit, and one of the main difficulties involved in the proof 
is to derive a semiclassical expansion with minimal regularity assumptions. The 
present article represents a summary of our recent work in [3] (see also [4j for a 
simplified model in one dimension) and we shall refer to [3| for all technical details. 

We emphasize that it is not rigorously understood how the BCS model approx- 
imates the underlying many-body quantum system. We will formulate the BCS 
model as a variational problem, but only heuristically discuss its relation to quan- 
tum mechanics. 

2. The Ginzburg-Landau Model 

The Ginzburg-Landau model is a phenomenological model for superconducting 
materials. Imagine a sample of such a material occupying a three-dimensional box 
If W denotes a scalar external potential, and A a magnetic vector potential, 
the Ginzburg-Landau functional is given as 



where ip € if 1 (A). Here, the Xi are real parameters, with A3 > 0. Note the factor 
2 in front of the A-field, which is reminiscent of the fact that if> describes pairs of 
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^For simplicity, we restrict our attention to three dimensions here, but a similar analysis applies 
in one and two dimensions as well. 
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particles. The microscopic model discussed in the next section does not have such 
a factor 2, as it describes single particles. 

By simple rescaling, we could take A3 = 2 1 A2 1 if A2 7^ 0. If A 2 > 0, one could then 
complete the square in the second line of dJ) and write it as Xs(l — \ip(x)\ 2 ) 2 — A3 
instead^ Since A2 can have either sign, however, we prefer to use the more general 
formulation in (fTJ) here. We shall derive formulas for the coefficients A^ from the 
BCS model below. 

The function ip is interpreted as the order parameter of the system. In the 
absence of external fields, i.e., for W — and A = 0, the minimum of (TTJ is 
attained at \tfj(x)\ 2 — ^/{^M) for A2 > 0, or at ijj = for A 2 < 0, respectively. 

Our main concern here will be the relation of (fT]) with the Bardeen-Cooper- 
Schrieffer theory of superconductivity, which we describe in the next section. We 
will derive GL from BCS in an appropriate limit, where the temperature is close 
to the critical one, and the external fields A and W are suitably small and slowly 
varying on a microscopic scale. 

We note that there is a considerable literature [6] concerning functionals of the 
type (fT]) and their minimizers, usually with an additional term added corresponding 
to the magnetic field energy. Such a term plays no role here since A is considered 
a fixed, external field. 



3. The BCS Energy Functional 

The BCS model is based on the microscopic many-particle Hamiltonian describ- 
ing the system under consideration. Consider a gas of spin 1/2 fermions confined 
to the box A. With /i denoting the chemical potential, the Hamiltonian is given by 

(2) H = J2((-^j + A( X] ) f - n + W( Xj )) +Y / V(\x. l -x J \). 

3 i<3 

Here, V is the two-particle interaction potential, which is taken to be local, i.e., 
it is a multiplication operator. In the original BCS Hamiltonian [2] it is replaced 
by an effective, non-local interaction which results from integrating out the phonon 
variables. For simplicity, we prefer to work with a local interaction of the form (|2|) 
hereH 

The BCS model can be obtained from @ by an approximation procedure (see 
Appendix A in [7]). One first restricts the states of the system to quasi- free states 
on Fock space. These states are completely determined by two quantities, the 
one-particle reduced density matrijfl 

(3) l{x,y) = (a f (x)a(y)) 
and the Cooper-pair wave function 

(4) a(x,y) = (a(x)a(y)) . 

One thus obtains an expression for the energy of such states solely in terms of 7 and 
a. One further ignores the direct and exchange term in this expression, and arrives 



2 This is, in fact, the convention used in [3]. 

^It is not necessary to assume V to be radial, as done here, but the formulas simplify somewhat 



in this case. 

4 We suppress the spin dependence in the notation for simplicity. 
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at the following BCS energy functional. With T > denoting the temperature of 
the system, 

■F(T) - Tr \((-iV + A(x)f- t i + W(x))j\ -TS(T) 

(5) + / V{\x-y\)\a(x,y)\ 2 dxdy. 

J AxA 

The BCS states can be expressed as 2 x 2 block matrices 

(6) r=f 1 " 

Y a 1 — 7 

satisfying the constraint < T < 1. The entropy 5(r) in ([5]) takes the usual form 
S(T) = -TiTlnT. 

One of the main questions concerning the model ([5]) is whether a minimizing state 
r has a identically zero or not, i.e., whether or not Cooper pairs exist. The non- 
vanishing of a implies a superconducting (or superfluid, depending on the context) 
behavior of the system. 

3.1. BCS Energy in the Translation-Invariant Case. It was shown in [7] that 
in the absence of external fields, i.e., for A = and W = 0, there exists a critical 
temperature T c > such that 

• T > T c : Minimizer is normal, i.e., a = and 7 = (l + exp((— V 2 — n)/T))~ l 

• T <T C : Minimizer has a^O. 

This critical temperature may be characterized by the operator 

(7) tf T „(-V*- rt + n W ), KrM^-^ 

having zero as the lowest eigenvalue. Note that the spectrum of Kt(—^7 2 — m) 
equals [2T, 00), hence is necessarily an isolated eigenvalue of (J7]) if T c > 0. 

This linear criterion on the critical temperature was used to derive precise asymp- 
totics of T c for weak coupling and/or low density, see [8]-[l0]. 

In the following, we shall assume that V is such that T c > 0, and that the eigen- 
function ao corresponding to the zero eigenvalue of ([7]) is unique. The potential V 
also has to be sufficiently regular to be form-bounded with respect to the Laplacian, 
and should decay at infinity. 

3.2. Microscopic vs. Macroscopic Scales. Let us introduce a small parameter 
h > 0, describing the ratio between the microscopic and macroscopic length scales. 
The external fields A and W occurring in the GL functional vary on the macroscopic 
scale, i.e., on the scale of the box A. The particle interaction V, on the other hand, 
varies on the microscopic scale. To take this into account, we replace the external 
fields A(x) and W(x) in the BCS functional © by 

(8) AO) = hA(hx) , W{x) = h 2 W{hx) 

and define the BCS functional on a rescaled box A = h~ 1 A. Here, x is the micro- 
scopic variable, while x = hx is the macroscopic variable. 

We find it more convenient to express the BCS functional in macroscopic vari- 
ables, and shall henceforth drop the ~'s. The resulting rescaled BCS functional 
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IS 

J"(r) = Tr + hA(x)f - /i + h 2 W(x)^j 7 -TS(T) 

(9) +/ VQi-^x-y^aix^dxdy, 

JAxA 

and it is defined on the /i-independent macroscopic volume A. Note the semiclassical 
nature of the appearance of h in the various terms in (|9]). For small h the energy 
is of order h~ 3 . 

In order for (J5|) to be well-defined, suitable boundary conditions have to be 
imposed on the boundary of A. In the limit h — > 0, these are not relevant in the 
bulk of the sample. In order to avoid technical problems related to these boundary 
condition, we chose in [3] to work with an infinite system instead, which is assumed 
to be periodic, and all energies are calculated per unit volume. Consequently, also 
the functions ip in the GL functional (JTJ) have to be periodic. 

4. Main Results 

We shall choose the temperature T in © to be close to the critical temperature 
defined in Subsection 13.11 More precisely, we take 

(10) T = T c (l - Dh 2 ) 

for some ^.-independent parameter D e R0 

Our main result concerns the asymptotic behavior of infr F(T) and the corre- 
sponding minimizers as h — > 0. 

Theorem 1. There exists a Ao > and parameters Ai, A2 and A3 in the GL 
functional such that 

(11) inf T(T) = F(T ) + X hM£ GL {ijj) + o(h) 

r 1/1 

as h — > 0, where Tq is the normal state (i.e., the minimizer in the absence ofV). 

Moreover, ifT is a state such that F(T) < J"(r ) + X h inf ^ £ Gh (ip) + o(h) then 
the corresponding Cooper pair wave function a satisfies 

(12) \\a - a GL \\h < o(h)\M\h = o(h)h- 1 
where 

(13) a GL (x,y) = h~ 2 ip y - ^ a (—J^j = Op(hifj (x)a Q (p)) 

and £ GL (^o) < m{^,S(tp) + o(l). 

Theorem[T]represents a rigorous derivation of Ginzburg-Landau theory. Starting 
from the BCS model, GL arises as an effective theory on the macroscopic scale, in 
the presence of weak and slowly varying external fields, and for temperatures close 
to the critical one. 

Remark 1. As mentioned above, -F(ro) is 0(h~ 3 ), hence the GL functional arises 
as an 0(h 4 ) correction to the main term. 



^The results in [3] were stated for D > 0, but their proof is equally valid for D < 0. 
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Remark 2. In (fl"3j) . Op denotes Weyl quantization. Theorem [T] demonstrates the 
role of the function ip in the GL model: It describes the center-of-mass motion of 
the Cooper pair wavefunction, which close to the critical temperature equals (|13[) 
to leading order in h. 

Remark 3. The method of our analysis can also be used to show that the GL 
model predicts the correct change in critical temperature in the BCS theory due 
to the external field [11]. As discussed in the next subsection, the parameter A2 
is proportional to the difference between the critical and the actual temperature 
(more precisely, to D in (fT0| ). and the mentioned shift corresponds to the largest 
A2 with the property that ip = is a GL minimizer. This can have either sign, 
depending on the external fields A and W. 

Remark 4. A similar analysis can be used at T = to study the low-density limit 
of the BCS model. In this limit, one obtains a Bose-Einstein condensate of fermion 
pairs, described by the Gross-Pitaevskii equation [T2l 1 1 3J . 

4.1. The Coefficients A^. The coefficients Ao, Ai, A2 and A3 in Theorem Q] can 
be explicitly calculated. They are all expressed in terms of the eigenfunction ao 
corresponding to the zero eigenvalue of ([7]), the critical temperature T c and the 
coefficient D in (fT0|). 

Specifically, if we denote by t the Fourier transform of 2Kx r ao, we have 

(14) A = J t{qf ( gi (f3 c ( q * fi)) + Wrfg^Ue ~ M))) ^3 > 

(15) A i = A o 1 i/ *(<z) 2 Si(/3c(9 2 -M)) ° ' 



D f x/_s2 -2 ( Pc, 2 d <l 



(16) A 2 = V W c L^ C ° Sh [W-tijW 
and 

(17) A 3 = A -^/ t{q) **W-»)) * 



16T,?y R 3 ^' g 2 -M (2tt)3- 

Here j3 c — T~ , and gi and gi denote the functions 

, , e 2z - 2ze z - 1 , . . 2e z (e z - 1) 
(!8) 9i(z)= z2(1 + e2)2 an d g 2 (z)= — - - 3 , 

respectively. One can show 3 that Ao > 0. Note that gi(z)/z > 0, hence also 
A 3 > 0. The coefficient A2 is proportional to D, defined in (ITD|) . 

As mentioned in Section [2J the terms in the second line of the GL functional 
(p} are often written as -^(1 — \ip(x)\ 2 ) 2 instead, with a suitable coupling constant 
K > 0. In our notation, n corresponds to 

(19) k=^/% 

(in case D > 0, i.e., T < T c ). 

Note that the normalization of ao is irrelevant. If we multiply ao by a factor 
A > 0, then Ao and A3 get multiplied by A 2 , while Ai and A2 stay the same. Hence 
the GL minimizer ipo gets multiplied by A -1 , leaving both \q£ gl (iJjq) and the 
product V'oao unchanged. In particular, ()13[) is independent of the normalization 
of ao- 
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5. Sketch of Proof 



In the following, we present a very brief sketch of the main ideas in the proof 
of Theorem Q] The actual proof is rather lengthy and we have to refer to [5] for 
details. 

The starting point is the identity 



(20) 

where T A = (1 + 
(21) H A = 



In 1 



,-Ha/t 



In 1 



,-Ho/T 



V(\x - y\/h)\a Gh (x,y)\ 2 dxdy 

T f 2 

-H(r, T A ) + / V(\x ~ y\/h) \a Gh (x, y) - a(x, y)\ dx dy 



,Ha/T\-1 



1 
A 



A_ 

-I 



f) = (-ihV + hA(x)) 2 -n+ h 2 W(a 



and A denotes the operator with integral kernel 

(22) A(x, y) = 2V(\x - y\/h)a GL (x, y) = 2hOp(Mx)(^V)(p)) ■ 
Moreover, H(T, Ta) denote the relative entropy 

(23) u(r, r A ) = Tr [r (in r - in r A ) + (i - r) (in (i - r) - in (i - r A ))] . 

It is non- negative, and vanishes only for T = T A . The latter property can be 
quantified in the form 



(24) 



T"tt(r,r A )>Tr K T (H A )(T-T A y 



The identity (|2U)) and the inequality (|2"4"|) hold true for any choice of the function 
ipo in (TT31). It is used in two separate steps. One first simply takes ipo = 0, and 
uses the gap of the operator Kt c (Ho) + V to conclude that a is close to an chgl 
for a suitable -0o, i.e., a is of the form (|13jl to leading order. One then repeats the 
argument with this choice of ipo, to conclude that the two terms in the last line of 
([20| are negligible for small h. 

To conclude the proof, it remains to calculate the terms in the first two lines 
in ([20)) . For this purpose one needs semiclassical estimates with good regularity 
bounds. The relevant estimates can be summarized as follows: 



Theorem 2. With errors controlled by H 1 and H 2 norms of ipo 
_ ^TTr [in (l + er H -/ T ) - In (l + e - H ^ T )' 
= h 2 V 2 (ip ) + h 4 V 4 (ip ) + \oh A £ GL (^ ) + 0(h 5 ) 



(25) 
and 

(26) 



V(\x - y\/h)\a Gh {x,y)\ 2 dxdy = h 2 V 2 {^ ) + h A V A {^o) + 0(h 5 ) 



for suitable T> 2 and T>±. 



In (|25|) we have already taken into account the choice ([10]) for the temperature T. 
For this choice, all the terms of order h 2 cancel, and one is left with Xoh 4 £ (ipo) 
after taking the difference between ([2"5|) and (|2"B"1) . 

This completes our very brief sketch of the proof of Theorem [T] 
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